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1. Introduction 

^f^ I Mankind has attempted to understand the origin of the Earth and planets 

cn ' for countless generations. In the past few decades, starting largely from the 

^^ [ pioneering work of Safronov in the 1960s, philosophical speculation about 

Qv ' their birth has been replaced with detailed analytical and numerical simu- 

Q> . lations whose goal is to reproduce the observed properties of our own solar 

(— I I system. With the recent groundbreaking discoveries of planets around or- 

Q^' dinary stars like our Sun (Mayor & Queloz 1995; Marcy &: Butler 1996; 

Butler & Marcy 1996; see also the review by Marcy in this volume), we 
^ ' now have a reasonable sample of nearby planetary systems that we can use 

c/3 . to test and validate theoretical models and from which we can learn new 

and unexpected features of planet formation. In the upcoming decade the 
launches of NASA's Space Interferometry Mission (SIM) and ESA's Global 
^ , Astrometric Interferometer for Astrophysics (GAIA) (and their longer term 

H I counterparts Terrestrial Planet Finder [TPF] and Infrared Space Interfer- 

ometry Mission [IRSI] ) will allow us to detect directly planets around nearby 
stars, which will give us an even more complete census of the number and 
content of planetary systems in the Solar neighborhood. To place this ob- 
servational data in a proper astrophysical context, we need to understand 
how planets are born and how they interact with their environment. In this 
review, I will discuss the dynamics of the growth of solid particles from 
initially micron-sized dust grains through kilometer-sized planetesimals to 
their final planetary mass. I will also discuss the formation of the gas-rich 
giant planets and their tidal interaction with the surrounding nebular disk. 
The reader interested in further details of the planet formation process 
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is encouraged to consult the excellent reviews by Lissauer (1993) and by 
Marcy and Butler (1998). 

1.1. OVERVIEW: FROM DUST TO PLANETS 

The low eccentricities and inclinations of the planets (see Table 1) have 
long been used as evidence that the Solar System formed from a flattened 
disk of gas and dust: the primordial solar nebula. The recognition that 
the collapse of magnetized, rotating molecular cloud cores naturally leads 
to protostellar disks with roughly the right mass and size (~ 100 AU) 
is one of the successes of modern star formation theory (for a thorough 
review see Shu, Adams, & Lizano 1987). The solid content in these disks 
is comprised of the interstellar grains that survive their passage into the 
disk (through a surface disk shock) plus nebular volatiles that condense as 
dust within the disk. In the current "standard model" of planet formation, 
the rocky cores of all planetary bodies grow from the accumulation of small 
dust grains; the increase in growth is tremendous - nearly a factor of 10^^ 
in size and lO^'^ in mass. In our Solar System, this model is sufficient to 
explain the formation of the terrestrial planets (Mercury, Venus, Earth, 
Mars), which have densities comparable to rock (see Table 1) and low-mass 
gas atmospheres that were most likely outgassed from the rocky material 
from which they accreted (Brown 1949; Prinn 1982; Pepin 1989). The giant 
planets (Jupiter, Saturn, Uranus, Neptune), which have much lower mean 
densities (c/. Table 1) and massive gas atmospheres, are also thought to 
have formed from initially rocky cores^. However, once the core mass grows 
to be larger than a critical value (^IOMq), rapid accretion of nebular gas 
transforms them into gas giant planets (Mizuno 1980). 

The growth of planet-sized solid bodies conventionally occurs in two 
loosely defined sequential phases. In the first phase (§3), dust grains, which 
are uniformly distributed in and strongly coupled to the nebular gas, grow 
through binary collisions. As they grow larger they settle gravitationally 
towards the nebular midplane forming a thin layer. Mutual collisions (and 
possibly gravitational instability of the dust layer) induce further growth in 
the midplane forming roughly kilometer-sized bodies called planetesimals. 
Planetesimals are, by definition, massive enough to have largely decoupled 
from the gas and move on nearly circular, Keplerian orbits about the cen- 
tral star. The first phase is believed to be relatively short, occurring on 
a timescale as little as a few thousand orbital periods (Weidenschilling & 
Cuzzi 1993). In the second phase, planetesimals continue to grow through 
inelastic collisions. The essential feature is that collision cross-sections are 

^Models invoking dynamical gravitational instabilities to form the giant planets 
(Cameron 1978) have fallen out of favor (see Lin & Papaloizou 1985). 
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TABLE 1 . Properties of Solar System Planets 
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Mercury 


0.39 


0.055 


5.4 


0.20 
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Venus 


0.72 


0.82 


5.2 


0.007 
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Earth (M®) 


1.0 


1.0 


5.5 


0.02 







Mars 


1.5 


0.11 


4.0 


0.09 
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Jupiter 


5.2 


318 


1.3 


0.05 
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Saturn 


9.6 


95.1 


0.70 


0.06 
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Uranus 


19 


14.5 


1.6 


0.05 




1 


Neptune 


30 


17.2 


1.6 


0.009 




2 



Notes: Pluto is not included because it is more typical of the low-mass 

outer Solar System Kuiper belt objects. 
IMfB = 6 X lO^^K = 3 X lO^^Mn 



increased above their geometrical value by the focusing nature of the gravi- 
tational interaction between the colliding planetesimals (Safronov 1972). If 
the relative velocity dispersion among the planetesimals remains sufficiently 
small, the increase in cross-section is greatest for the most massive plan- 
etesimal. This planetesimal will then collide more frequently and gain mass 
more rapidly than its neighbors, which causes its collision cross-section to 
increase even further. "Runaway" growth occurs in which the most massive 
planetesimal accretes all the available solid matter in its local "accretion 
zone" - the region of the disk over which its gravity can perturb nearby 
planetesimals into colliding orbits. The runaway process produces a series 
of rocky protoplanetary bodies that are dynamically and physically isolated 
from each other and occurs on timescales of ~ 10^ years (Wetherill & Stew- 
art 1989; Aarseth, Lin, & Palmer 1993; Weidenschilhng, et al. 1997). After 
the runaway stalls because all the available nearby solid matter has been ac- 
creted, further collisional growth of the protoplanets can occur but only on 
much longer timescales (~ 10^ years) as mutual gravitational perturbations 
cause orbits to cross and highly inelastic collisions to occur (Wetherill 1980; 
Chambers, Wetherill, & Boss 1996). 

As the planets grow, they find themselves in a nebular environment 
that changes with time as the gas in the disk is accreted by the central star 
(Ruden &: Lin 1986). Eventually, the disk is cleared of gas via processes 
such as stellar winds or UV irradiation by the central star (Shu, Johnstone, 
& Hollenbach 1993). The growth of the terrestrial planets to their cur- 
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rent masses could have continued long after the dispersal of the gas disk. 
However, because the more massive atmospheres of the giant planets were 
accreted from the nebular gas disk, the growth of their rocky cores to the 
critical mass had to be complete before disk dispersal. Observations indicate 
protostellar disk lifetimes are ^10^ years (Strom, et al. 1989; Beckwith, et 
al. 1990; Zuckerman, Forveille, & Kastner 1995; see also the review by Beck- 
with in this volume), placing an important time constraint on the planet 
formation process. Current models have difficulty forming the giant planet 
cores on this timescale unless the mass of the solar nebula is significantly 
larger than the minimum-mass solar nebula. 

The giant planets are massive enough that gravitational tidal interac- 
tions between the planet and the surrounding disk can create annular gaps 
largely devoid of gas at the orbit of the planet. This mechanism of tidal trun- 
cation (Lin & Papaloizou 1993) has been invoked to explain why the giant 
planets stop their rapid gas accretion and thus attain their final masses 
(§5.2). The tidal interaction also causes the orbital semimajor axes of the 
planets to decrease in time. This orbital migration of the giant planets is a 
serious problem that will be discussed in more detail in §5.3. 

1.2. FUTURE RESEARCH QUESTIONS 

Throughout this review I will list a selection of the major unsolved ques- 
tions in planet formation theory. I hope these questions will provoke and 
encourage the reader to consider devoting some thought to their solution. 

— When does the planet building epoch begin? 

• Is it a continuous process beginning from the initial formation of 
the disk? 

• Does planet formation require some critical conditions (such as a 
sufficiently cool or quiescent nebula) before it begins? 

— How does planet formation modify the disk environment and its evo- 
lution? 

• Are their unambiguous observable signatures such as holes or 
gaps? 

— What outcomes are likely? 

• A diverse planetary system like our own? 

• Giant planets close to the central star as in 51 Pegasus? 

• Rubble disks like (3 Pictoris? 
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2. Structure and Evolution of Disks 

Protoplanetary disks are just one example of an accretion disk - a geomet- 
rically thin disk of gas (and dust) revolving around a central mass point 
(Pringle 1981; Lin & Papaloizou 1985; and the review chapter by Kenyon 
in this volume). We can understand their basic properties by examining 
the force balance equations assuming the gas is in orbit around a central 
protostar with mass M* . Adopting cylindrical coordinates (r, (p, z) and as- 
suming axisymmetry, the disk is in vertical hydrostatic equilibrium where 
the vertical pressure force balances the vertical component of the central 
star's gravity 

19P CAWA^_„.^, (1) 



p dz r^ 

where P and p are the gas pressure and density, and where we have defined 
the Keplerian rotation rate 



We can scale equation (1) to find an expression for the vertical height of 
the nebula, H, 

where the nebular sound speed c is defined by the relation P = p(?. We 
say the disk is geometrically thin \i H <^ r, which from equation (3) is 
equivalent to the disk being dynamically cold, c <C f^r, i.e., the nebular 
sound speed is much less than the orbital Kepler speed. Because disks are 
thin, it is convenient to use vertically averaged quantities such as the surface 
density, S: 

1^ = £" pdz^2pH = 2^, (4) 

where we interpret p in the last two terms on the right side of the equation 
as the value of the density at the midplane of the nebula. 
Force balance in the radial direction is given by 

v^ GM^ IdP 
r r^ p or 

where v is the orbital speed of the gas. For thin, cold disks, the radial pres- 
sure gradient is smaller than the stellar gravitational force by ~ (H/r)'^ <^ 1 
(as may be verified by scaling eq. [5] and using eq. [3]) and is typically ne- 
glected in most accretion disk applications. In this case, equation (5) shows 
the gas rotates in centrifugal balance at the Keplerian orbital speed, v = Or. 
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In most disks, the pressure is largest near the star and decreases outward 
so that the radial pressure term in equation (5) is negative. This means the 
pressure of the gas partially supports the disk against gravity, and the disk 
gas rotates at a slightly sub-Keplerian speed v = i}r — Av with At; > 0. If 
we substitute this expression for v into equation (5), we find 



1 

At; = - 
2 



dlnP 



dlnr 






where we have used equation (3) in the last terms on the right. The depar- 
ture from Keplerian motion is indeed very small being much less than the 
sound speed which is itself much less than the Keplerian rotation speed. 
We shall see in §3.1 that this small departure from exact Keplerian rotation 
plays a critical role in the motion of solids in protoplanetary disks. 

In two seminal papers, Lynden-Bell &: Pringle (1974) and Shakura & 
Sunyaev (1973) demonstrated that the evolutionary behavior of accretion 
disks is governed by the outward transport of angular momentum, which 
allows mass to flow inward to be accreted by the central object. The precise 
mechanisms that cause the transport in any astrophysical disk are still far 
from certain. Numerous suggestions have been made in the case of protostel- 
lar disks (see the review by Adams &: Lin 1993): gravitational instabilities 
(Lin & Pringle 1987; Shu, Tremaine, Adams, & Ruden 1990), magnetic 
instabilities (Stepinski &: Levy 1990; Balbus & Hawley 1991), thermal con- 
vective instabilities (Lin & Papaloizou 1980; Ruden & Lin 1986), stellar 
or disk winds (Hartmann & MacGregor 1982; Wardle & Konigl 1993) and 
fluid dynamical shear instabilities (Dubrulle 1993). In most applications, 
the details of the physical cause of the transport are ignored and angular 
momentum is assumed to be transported by some form of localized fluid 
dynamical turbulence that creates an eddy viscosity in the gas. The disk 
evolution is governed by the magnitude of the turbulent viscosity, v, which 
is the product of the turbulent eddy size, it, and velocity, Vt- The eddies are 
assumed to move subsonically and to have sizes smaller than the vertical 
height of the disk, which leads to the well-known parameterization of the 
viscosity as (Shakura & Sunyaev 1973) 

v rsj y^ . i^ = acH, (7) 

where a;$l is a dimensionless parameter. Under the action of viscous stress, 
an initial surface density distribution evolves diffusively by spreading out- 
ward while transferring mass inward to be accreted by the central protostar. 
The viscous evolutionary timescale is 

r^ 1 1 /r \2 
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. 10' C^] m)' (^^r- years. 

The sohds in protoplanetary disks find themselves in a continuously 
changing nebular environment in which temperature, density, pressure, and 
turbulent velocity typically decrease with time (after transients from the 
initial distribution of matter decay; see Fig. 1). Detailed calculations using 
reasonable values for the eddy viscosity (a ~ 10~^) yield evolutionary times 
of order 10^ - lO'^ years (Ruden & Lin 1986; Ruden k Pollack 1991), which 
are in agreement with the observational disk lifetimes cited above. 

2.1. THE MINIMUM-MASS SOLAR NEBULA (MMSN) 

The Solar System was formed from the collapse of a molecular cloud having 
solar element abundances, yet the present day planets are clearly of non- 
solar composition. We can ask the important cosmogonical question, "How 
much solar composition gas is required to make the planets?" The mass 
and inferred spatial distribution of this gas are known as the Minimum- 
Mass Solar Nebula (hereafter, MMSN). The answer to this question has 
been given by several authors (Weidenschilling 1977a; Hayashi 1981) who 
take the present amount of condensed solid matter in the Solar System 
and augment it to solar composition to determine the minimum mass. The 
radial surface density distribution of this mass is deduced by smearing 
out the augmented mass of a planet over an area defined by its nearest 
neighbors.^ 

A key ingredient in the analysis is the dust-to-gas ratio, Ci which mea- 
sures the mass fraction of dust (i.e., high-Z elements) in a solar compo- 
sition gas. In cool regions of the nebula with temperatures below about 
170 K, both rocky and icy material (such as water, methane and ammonia) 
can condense from the gas giving the dust-to- gas ratio a value C ~ 1/60 
(Hayashi 1981). In warmer regions (above 170 K but below about 1500 K 
where much of the dust evaporates) only the more refractory (high con- 
densation temperature) solids survive and the dust-to-gas ratio drops by a 
factor of four to (^ ~ 1/240 (Hayashi 1981). To find a rough value for the 
mass in the MMSN we must estimate the present amount of solids in the 
planets. This comes almost entirely from the high-Z material inside the 
giant planets, which is estimated to range from 40 - 80 M® (Hubbard & 
Marley 1989; Chabrier, et al. 1992) compared to only 2 M® in the terres- 
trial planets. Adopting a value of 60 Mq for the present mass of condensed 
solids and an average dust-gas-ratio of 1/100 gives a rough estimate of 6000 
M® ~ 0.02Mq for the mass of solar composition gas needed to make the 

^An important assumption is that the planets formed at their present locations. 
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Figure 1. Models for the evolution of the primordial solar nebula. Surface density, E, 
aspect ratio, H/r, central (midplane) temperature, Tc, and surface (effective) tempera- 
ture, Ts, are plotted versus distance from the sun (in AU) for the times a) 10"', b) 10**, 
c) 10^, d) 10^ years. Note the diffusive spreading of the nebula from less than 10 AU to 
100 AU, and the cooling of the nebula with time. The disk always remains geometrically 
thin {H/r <C 1), and at late times H/r « 1/25. The surface temperature falls off as r~ ' 
(marked as the dotted line) away from the outer disk edge. The viscous a parameter is 
10^^ in these calculations adapted from Ruden & Lin (1986). 



planets.^ A more detailed calculation of the MMSN by Hayashi (1981) is 
illustrated in Figure 2. Disk masses derived from millimeter observations of 
nearby star forming regions are found to range from 0.005 - 0.2 Mq (Beck- 



^Adding the mass of solar composition gas already present in all the planetary atmo- 
spheres (;S4OOM0) is only a small correction. 
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with & Sargent 1996; also see the review by Beckwith in this volume), 
comparable to the MMSN. 
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Figure 2. Surface density distribution in the Minimum-Mass Solar Nebula (MMSN) 
adapted from Hayashi (1981). The planets are denoted by Me, V, E, Ma, J, S, U, N. 
Outside 2.7 AU the temperature is deduced to be cool enough for ices to condense, which 
causes the surface density of solids to jump by a factor of four. The gas distribution 
shows no discontinuity, and the surface density of both gas and solids falls off as r~^'^. 
The total mass between 0.35 AU and 36 AU is 0.013 M©. 



The MMSN, while playing a useful role as a conceptual guide, is perhaps 
becoming overused as a quantitative "reference" standard in model calcu- 
lations of the planet formation epoch. The planets formed in an evolving 
disk whose properties (such as the position where ices condense) changed 
with time; the static r"^'^ surface density distribution presented by the 
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MMSN may not only be quantitatively wrong, but misleading is well. If it 
does accurately describes any part of the evolution of the solar nebula, it is 
probably the final stages when planet formation is complete rather than the 
initial stages when planet formation began. Furthermore, as we will see in 
§5.3, it is now thought that orbital evolution of the giant planet semimajor 
axes may have occurred, which would modify any current reconstruction of 
the MMSN surface density distribution from that depicted in Figure 2. 

2.2. NONLINEAR NATURE OF DISK PHYSICS 

It is obvious there is a very complicated interrelationship among all the 
physical processes that can occur in protoplanetary disks. Figure 3 is an in- 
tentionally complicated schematic diagram indicating some of the predom- 
inant physical processes, their dependencies and their feedbacks. Although 
there are many models by many different research groups addressing the 
key underlying physics in any one of the boxes in Figure 3, successful inte- 
gration of the different processes together into a coherent whole has been 
elusive. There are several reasons for this: the inherent complexity of the 
physics involved, the lack of experimental guidance as to the magnitude of 
important physical parameters, the inability of current observational tech- 
niques to probe the relevant scales (which occur at angular resolutions of 
~ 10"^ arcseconds in the nearest star forming complexes) and the nonlinear 
interplay among the processes themselves. For these reasons, we are far from 
having a theory that can predict the type of planetary system formed from 
a given nebular environment. Furthermore, researchers have typically fo- 
cused their modeling on explaining the properties of our own Solar System, 
perhaps rejecting scenarios that, although inadequate for our own system, 
might be good representations of the physics in other protoplanetary disks. 
As is often the case, more complete observational data will enable modelers 
to explain the diversity of planetary systems found in nature. 

3. From Grains to Planetesimals 

3.1. PARTICLE-GAS DYNAMICS 

The motion of solid particles in protoplanetary disks is affected by the 
drag they experience as they move through the gaseous background of 
the nebula (Whipple 1972; Adachi, Hayashi, & Nakagawa 1976; Weiden- 
schilling 1977b). We will consider the solids as spheres of radius a having 
internal mass density 5 (to be distinguished from the ambient gas density 
p) and mass m = AirSa^ /3. The form of the gas drag force, Fd, depends on 
the size of the particle, its speed through the gas, and local gas properties. 
The Epstein drag law, which is appropriate for particles with sizes smaller 
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RELATIONSHIPS AND DEPENDENCIES 
AMONG PROTOPLANETARY DISK PROCESSES 

Figure 3. Schematic diagram showing the relationships and feedbacks among nebular 
processes, many of which are only poorly understood. Small modifications to the under- 
lying physics for any one process may have profound effects on the ultimate evolutionary 
outcome of the star/disk system. 



than the mean free path in the gas, is 



Fd 



--vra /)cv, 



(9) 



where c is the gas sound speed and v is the velocity of the particle relative 
to the gas. In most regions of the solar nebula, this law is valid for particles 
with sizes smaller than a few centimeters (see the references above for a 
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complete list of the different drag law regimes) ; for simplicity of presentation 
in this review, I will adopt equation (9) as the principal drag formula. The 
stopping time for a particle is defined to be the time needed for drag to 
dissipate its momentum, tg = mv/Fi). Substituting equation (9) into this 
definition and eliminating the volume gas density in favor of the surface 
density from (4), we find the relation 

Ots = ^. (10) 

The dimensionless quantity 0,ts is the ratio of the particle stopping time to 
the local orbital time in the disk and is a key determinant of the dynamical 
behavior of the particle. In the limit 0,ts <C 1, the particle, considered 
"small", is very strongly coupled to the gas. Drag will force it to move 
at the local gas velocity, which we saw in §2 to be sub-Keplerian. In the 
opposite limit, fits 3> 1, the particle, considered "large", is only weakly 
perturbed by gas drag and orbits at the Keplerian speed. 

Let us consider the dynamics in more detail. At each radius, gas drag 
forces small particles to orbit at the local gas speed, which is smaller than 
Keplerian. The net effective gravity (stellar plus centrifugal) acting on the 
particle points toward the central star, and the particle drifts inward at the 
terminal speed. This radial velocity is easily shown to be 

Vr = -2ntsAv, (17ts < 1) (11) 

where Av is the difference between the Kepler and gas speeds defined in 
equation (6). Large particles move at the Kepler velocity and experience a 
headwind as they move through the more slowly moving gas. The gas drag 
torque causes the orbit of the particle to decay with radial velocity 

2Au 
Vr = —^- (^ts > 1) (12) 

Detailed calculations by Weidenschilling (1977b; see Fig. 4) verify these 
limiting behaviors and show that the maximum inward drift speed is ~ At; 
and is achieved by particles with sizes that satisfy fits ~ 1 • There are several 
points to note about this inward drift. First, it is size dependent leading to 
differential speeds among particles with different radius, which can enhance 
the rate of collisions. Also, the relative speed between bodies entrained 
in the flow can easily exceed their escape velocities (see Fig. 4). Second, 
the affect is largest for particles having radii for which 0,ts ~ 1. Using 
equation (10), we find that the size of these particles is a ~ S/5, which is 
of order 1 m at the position of the Earth. Third, although the magnitude 
of the drift speed is highly subsonic, particles satisfying Ots ^ 1 decay into 
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the central star in a time very much less than the evolutionary time of the 
disk. The drift time for these particles is tdrift = r/vr = ^~^ {'''/H) ~ 
atvis ~ 100(r/AU)^/^ years. Of course such particles are likely to collide 
with other solids on their way in and through shattering or growth move to 
size regimes with smaller drift rates. Nonetheless, gas drag-induced decay 
is a potentially large sink for the solid component of the disk. 
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Figure 4- Radial and transverse velocity components, relative to tlie gas, as a function 
of particle radius. The difference between tlie Kepler speed and the gas speed is labeled 
AV. Small particles have transverse speeds equal to the gas speed, while large particles 
move at the Kepler speed, i.e., have transverse speeds relative to the gas of AV. The 
maximum inward drift speed is « AF and is achieved by particles with Qts ~ 1. The 
changes in slope are transitions between drag regimes, and Va is the escape speed from 
the particle. From Weidenschilling (1977b). 



3.2. MIDPLANE SETTLING AND GROWTH 



As discussed earlier, the formation of solid bodies in the midplane of the 
nebula begins with agglomeration and sedimentation of smaller grains ini- 
tially suspended throughout the nebula. Small particles settle in the verti- 
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cal gravitational field {gz = {z/r)GM^/r'^ = il^z) at the terminal velocity. 
Using the Epstein law (9), the settling speed and the time to reach the 
midplane are (fitg ^ 1) 

vz = -{nts) nz, (13) 

The last equation shows that the settling time is always much longer than 
the rotation period for small grains (Otg ^ 1) and that small grains set- 
tle more slowly than larger grains (tgcttic oc 1/a). For typical solar neb- 
ula conditions the settling time for micron-sized grains is about 10^ years, 
comparable to the lifetimes of protoplanetary disks. Without collisions to 
increase grain size (leading to more rapid settling), small grains, which are 
the dominant contributor to the opacity, will remain suspended in the disk, 
which will remain optically thick. 

Safronov (1972) was one of the first to realize that collisions not only 
make particles grow but also rapidly decrease the settling time. Equa- 
tions (13) and (10) show that the settling speed is proportional to the 
size of the particle. Larger grains will sink towards the midplane faster and 
will sweep up smaller grains in their path. This causes them to grow larger 
and to sink even more rapidly. In this manner, the dust component of the 
disk will therefore "rainout" and form a thin layer of much larger bodies in 
the midplane of the nebula. We can estimate the size of these bodies by us- 
ing mass conservation. The total surface density of solids in the disk is C'^, 
where C, is the dust-to-gas ratio introduced above. A body that reaches the 
midplane with final radius a will have swept up a column of dust having a 
mass ps • CS • 7r(a/2)^, where we have taken its average radius on the descent 
to be a/2 and have introduced the probability of sticking in a collision, pg. 
Equating this to its mass m = 47r(5a^/3 yields a size 

which is about 1 cm at the Earth in the primordial solar nebula (for ps = 
1). The time for solids to sediment to the midplane can be estimated by 
inserting this value for a into equation (14); more accurate integrations 
yield (Nakagawa & Hayashi 1981) 

_ 100 J_ 

'■sediment ^~^ ~pz~ 7) \^^) 

n psC 

which is only a few thousand orbital periods at any point in the disk (again 
for Ps = 1). 
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There have been several detailed numerical simulations of the collisional 
growth and sedimentation of dust in the solar nebula (Weidenschilling 1984; 
Nakagawa, Nakazawa, & Hayashi 1981; Mizuno 1989; Weidenschilling & 
Cuzzi 1993). A major difficulty with simulating particle growth is our in- 
adequate knowledge of the physical inputs. The sticking probability is a 
complicated and currently only poorly understood function that depends 
on the impact velocity and the size, shape and internal strength of the 
impactors, but it is usually simply taken to be unity for all collisions in 
the simulation. Realistic nebular particles may have low internal strengths 
and may shatter in high speed impacts. Some experimental work has been 
performed on this important issue. Blum & Miinch (1993) did not find any 
sticking in collisions between millimeter-sized dust aggregates at speeds 
^1 m/s and found fragmentation at speeds above this. More recently, this 
group (Wurm & Blum 1998) has investigated lower speed (^1 cm/s) colli- 
sions between small, porous, fractal grains with sizes ~ 10 /um and found 
perfect sticking (ps = 1). In the outer, cooler regions of the nebula volatile 
ices can form frosts on particle surfaces, which may enhance the sticking 
probability (Bridges, et al. 1996; Supulver, et al. 1997). 

Turbulence in the gas can have profound affects on the evolutionary 
simulation. A given sized particle will respond to the fluctuating turbulent 
velocity field according to whether its stopping time is longer or shorter 
than the local eddy turnover time. Stirring by turbulence can inhibit sedi- 
mentation of particles that have sizes small enough that they are strongly 
coupled to and co-move with the turbulent eddies. The turbulent velocity 
field can also greatly increase the relative velocity between particles with 
different sizes (Markiewicz, Mizuno, & Volk 1991), which can increase the 
collision rate but also lead to collisions where the bodies shatter (Weiden- 
schilling 1984). Given our poor understanding of turbulence in protoplan- 
etary disks, its effects on grain growth (and vice versa) and grain sedi- 
mentation are presently far from clear. Even with these uncertainties and 
although the individual details vary, the different numerical simulations 
nonetheless yield midplane particle sizes and sedimentation timescales that 
are not too dissimilar from the simple estimates above (eqs. [15] & [16]). 

3.2.1. Radial Distribution of Solids in Protoplanetary Disks 

The sedimentation of dust grains produces a thin layer of solids in the 
nebular midplane. The surface density of this distribution is often taken 
to be just the gas-to-dust ratio times the local gas surface density, (T,. To 
test whether this is a valid approximation and to see whether the solid 
surface density distribution of the MMSN is reproduced, I have calculated 
simplified models of the midplane accumulation of solids using the nebular 
evolutionary models of Ruden & Pollack (1991). The models follow track 
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the condensed and vapor phases of refractory rocky and volatile icy mat- 
ter, with ices condensing only in nebular regions with temperatures below 
170 K. Early in the evolution, the nebula is too hot for ices to condense 
anywhere, although rocky matter can condense everywhere outside about 
0.1 AU. As the nebula evolves and cools, the radius where ices can con- 
dense, i.e., the T = 170 K boundary, moves inward and icy solids sediment. 
The treatment of grain growth and sedimentation is phenomenologically 
modeled by allowing the condensed solids to rainout on a timescale given 
by tq/CQ {c.f. eq. [16]), where tq is an adjustable constant. This model 
is nonlinear because as the solids sediment to the midplane, the dust-to- 
gas ratio C decreases, and the sedimentation time lengthens. Figure 5 shows 
the results of two evolutionary scenarios: the first with rapid sedimentation. 
To = 600, and the second with much slower sedimentation, tq = 6 x 10'^. 
Both models started from the same initial conditions. Significantly more 
solids are deposited in the midplane in the rapid scenario, but at the time 
shown in the plot (10^ years) the nebula is warm enough that ices only con- 
dense outside 6 AU. The results of the slow scenario give a surface density 
distribution comparable to the MMSN after ~ 10® years, but with a more 
shallow slope (and more solid matter) in the outer cool regions. 

It seems clear from these idealized models that very different distri- 
butions of midplane solid matter are possible, from the same initial mass 
reservoir, if the sedimentation physics (primarily the average sticking prob- 
ability) is different. The solid distributions shown in Figure 5 do not bear 
a simple constant multiplicative relationship (C^) to the instantaneous gas 
surface density distribution. The reason for this is the solids at any given 
radius were deposited over time and reflect accumulation from different neb- 
ular conditions, with primarily rocky material sedimenting early on when 
the temperature was high followed by icy sedimentation when the nebula 
was cooler. Note also that the locus at which ice can condense is different 
in the two models; in general, the midplane distance at which T = 170 K 
(the ice condensation front) is a function of time. 

3.3. PLANETESIMAL FORMATION 

Planetesimals are bodies large enough to move on Keplerian orbits largely 
unaffected by gas drag forces over the age of the disk. For typical nebular 
conditions, bodies larger than about 1 km satisfy this constraint. A key 
unanswered question is how the smaller sized bodies (~ 1 cm; eq. [15]) pro- 
duced in the sedimentation process assemble into planetesimals. The most 
elegant and efficient mechanism proposed is that the thin layer of solids pro- 
duced by sedimentation becomes gravitationally unstable and rapidly forms 
fragments that collapse to produce planetesimals (Goldreich & Ward 1973). 
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Figure 5. Radial distribution of solids in a model of grain growth and sedimentation 
(Ruden 1998). Surface density (in g/cm^) is plotted versus distance in AU. The surface 
densities of icy (dotted), rocky (dashed), and rock + ice (solid) matter are plotted. The 
MMSN surface density distribution is shown for comparison (dot-dashed). Panel (a) [(b)] 
is for the rapid [slow] accumulation scenario, with the evolution times indicated. Note 
the difference in vertical scale between the panels. See text for further details. 



They showed the critical fragmentation lengthscale, Ac, is 



\r ~ A-irr 



M, 



and the typical fragment mass, Mqw is 



Mgw = CSA^ = levrM, 






(17) 
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In these equations, it is the surface density of solids, CS, that is important, 
and the ratio in parenthesis is roughly the mass of the sohd disk to the 
protostehar mass. Using nebular parameters from the MMSN, Goldreich & 
Ward (1973) showed that kilometer-sized bodies can be produced by this 
gravitational collapse process in a few thousand years. 

The Goldreich- Ward instability occurs only in very quiescent dust disks. 
If the random velocities of the particles composing the dust disk are too 
large (or equivalently if the dust layer is too thick)^ then the instability 
is quenched. The critical random speed is vrGCS/il, which is only about 
10 cm/s. In the thin dust disk the solids rotate at nearly the Keplerian ve- 
locity while the gas above the layer rotates at a sub-Keplerian rate (eq. [6] ) . 
It is thought that turbulence generated by this shear layer will stir up the 
random velocity in the dust disk and prevent the Goldreich- Ward instabil- 
ity (Weidenschilling & Cuzzi 1993). Detailed hydrodynamical models of the 
particle-gas interaction have been made by Cuzzi, Dobrovolskis, & Champ- 
ney (1993) who find that it is unlikely the dust layer can achieve such low 
random speeds. 

If the Goldreich- Ward mechanism does not proceed in protostellar disks, 
then it is likely that planetesimals are produced in the thin dust layer by 
inelastic binary collisions. Recall that gas drag causes radially inward flow 
of solid particles (eqs. [11] &: [12]). In the dense particle layer, solids domi- 
nate the local density and the drift speed is reduced (Cuzzi, Dobrovolskis, 
& Champney 1993). Since larger particles have lower drift speeds (eq. [12]), 
Cuzzi, Dobrovolskis, & Champney (1993) have argued that planetesimals 
can grow in this layer as the smaller solids drift past them and are collision- 
ally accreted. They estimate that growth to 100 km sizes can be achieved 
in ^10^ years. 

3.4. PROSPECTS FOR OBSERVING THE EARLY PLANET FORMATION 
EPOCH 

It would be satisfying if theoretical models gave an unambiguous prediction 
for the observational signature of grain growth in young protoplanetary 
disks, for example, the existence of a "hole" in the infrared spectrum (such 
as the one inferred observationally in HR 4796 by Koerner, et al. 1998). 
Unfortunately, this is presently not the case. In order to "see" forming 
planets or other structures in the midplane of the disk, the optical depth 
to the nebular midplane must be small. The disk optical depth is equal 
to the gas surface mass density times the opacity. The midplane regions 
can remain cloaked to the eyes of observers if either is sufficiently large. 

^The factor in parentheses in equations (17) & (18) is the maximum value of the 
aspect ratio, H/r, the layer can have to be gravitationally unstable. 
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The primary contributors to the opacity are small grains, but its value 
even for the standard interstellar grain-size distribution is subject to debate 
(Beckwith &: Sargent 1993; see the review chapter by Beckwith); grain size 
evolution only complicates the picture (Miyake & Nakagawa 1993). 

As noted in §3.2, micron-sized grains take times comparable to the age 
of disks to settle in quiescent nebulae. Any small grains left suspended in 
the disk after the bulk of the solid mass has sedimented to the midplane 
will remain suspended in the disk as long as it is quiescent. The effects of 
nebular turbulence on a suspended small grain population are harder to 
discern. By itself, turbulence causes small particles to remain stirred up 
and suspended in the nebula. However, turbulence also mixes these parti- 
cles into the midplane region where they may collide and stick to larger 
bodies (Weidenschilling & Cuzzi 1993). The net result of these two com- 
peting processes on the small grain population is unknown. Also uncertain 
is how the turbulent properties change as the disk optical depth changes, 
i.e. whether the disk will remain turbulent as grain growth proceeds. 

Current theoretical models cannot reliably predict how much of the ini- 
tial solid mass will remain leftover as small particles, from processes such as 
particle shattering and erosion. The nebular models of Ruden & Lin (1986, 
Fig. 1) and Ruden & Pollack (1991) remain optically thick through mid- 
infrared wavelengths at all disk radii as long as ^1 % of the initial grain 
mass remains suspended. The MMSN (Fig. 2) will also remain optically 
thick out to ^10 AU under the same conditions. Finding observable signa- 
tures may require waiting until the disk gas is cleared (probably taking with 
it the small grain population) or searching for planet formation signatures 
in long wavelength regimes that remain optically thin. 

3.5. FUTURE RESEARCH QUESTIONS 

— What are the appropriate conditions in protostellar disks at the incep- 
tion of the planet building epoch? 

• Do they resemble the MMSN? 

• Do different initial conditions lead to diverse planetary systems? 

— What is the sticking probability in grain-grain collisions? 

• Do simulations using Ps = 1 adequately model grain growth? 

— Is planet formation inefficient? 

• How much solid mass is lost via gas-drag induced radial flow? 

— How are planetesimals formed? 

• Via the Goldreich-Ward instability or by collisions? 

— What are observable signatures of planet formation? 
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• Holes, or gaps, or warps? 

4. From Planetesimal to Planet 

Planets are formed from the midplane planetesimal distribution by inelas- 
tic collisions. The key feature of the collisional evolution is that the cross- 
section for pairwise collisions can be increased above the geometrical value 
by gravitational attraction between the colliders (Safronov 1972). This re- 
sult can be shown quite easily by considering energy and angular momen- 
tum conservation in the collision between two particles with masses mi and 
1712 and radii ai and 02. In the center-of-mass system, the constant orbital 
energy and angular momentum are 

1 2 Gmim2 1 2 f-,n\ 

= 2^"^ r "2^ ' ^ ^ 

J = ^ |r X v| = fj,rvsm6 = fj,bV, (20) 

where /U = mim2/{mi + 171,2) is the reduced mass. In the rightmost equali- 
ties I have evaluated the constants E and J when the particles are at large 
separation where b is the impact parameter and V is the relative collision 
velocity. Without gravity the maximum impact parameter that results in a 
collision is the sum of the radii ai + a2- With gravity, a grazing collision re- 
sults when the minimum separation isr = ai + a2 and the velocity is purely 
tangential, 6 = 7r/2. Substituting these two conditions into equation (20) 
gives the speed at impact as v = bV/{ai + 02), which can be substituted 
into equation (19) to give the gravitationally enhanced cross-section for 
collision 



7r6^ = TT (ai + 02) 



where the mutual escape speed, Ve, is 



1 H 



(21) 



.e = J'-^iH^ll^. (22) 

V ai + 02 

The term in brackets in equation (21) is the gravitational focusing factor 
and can greatly increase the cross-section above geometrical if the relative 
collision speed, V, is much smaller than the escape speed, Vq. Also note that 
if the focusing factor is large the cross-section is proportional to the fourth 
power of the particle radius (because v"^ oc a^) rather than the second. 
In honor of Safronov's pioneering contributions to the study of collisional 
planet growth, the gravitational focusing factor is often written as 1 -|- 29, 
where 9 = v1/2V'^, is called the Safronov number. Detailed three-body 
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numerical orbit integrations confirm that the cross-section (21) is a good 
approximation (Ida & Nakazawa 1989; Greenzweig & Lissauer 1990, 1992). 
The equation for the collisional growth of a planet with mass rrip = 
47r(5a^/3 (recall 6 is the internal density of the planet) accreting matter 
from a background "swarm" of planetesimals that move with respect to it 
at a relative speed V is (using the cross-section [21] and neglecting the size 
of the planetesimals) 



dnip 2 

—-f- = ps-V -TTa 
at 



l-\ — - 



(23) 



where ps is the volume density of the planetesimal swarm. We can rewrite 
this as an equation for the rate of change of the planetary radius 



da 1 /ilS, 



dt 8 \ 6 



1 H 



(24) 



where Sg = Ip^ jQ. is the surface density of the planetesimal swarm (c./. 
eq. [4]). 

4.1. ORDERLY GROWTH 

The parameter that plays the most crucial role in determining the type of 
planetary growth is the Safronov number, i.e., the magnitude of the plane- 
tary escape speed relative to the planetesimal velocity dispersion. Using a 
variety of analytical approaches, Safronov (1972) argued that as planetes- 
imals accumulated, gravitational scattering off the largest member of the 
swarm would cause the velocity dispersion to keep pace with the escape 
speed from that member, regulating Q to be near unity. In this scenario 
planets grow in an orderly fashion with most of the mass in the large bod- 
ies. From equation (24), the accumulation time is 

a 1 f 6 a\ ,^ _ ,„„, 

*-'*"' = *^ " n ( et) • <* ~ '> P^' 

" i«''(^)%-- 

The growth time is determined by both the Kepler clock and the solid 
surface density. In the last line I have evaluated the growth time for Earth- 
sized bodies and taken Eg from the MMSN distribution (Fig. 2). These 
timescales are quite long, far greater than the lifetime of nebular disks, 
with the orderly growth process being unable to produce the cores of the 
outer giant planets in less than the age of the universe! 
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4.2. RUNAWAY GROWTH 

The key to forming planets more rapidly is keeping the velocity dispersion 
small enough so that the gravitational focusing factor becomes very large 
{9 ^ 1). In this case, the most massive planetesimal in the swarm will have 
the largest collisional cross section (oc a^), will grow more rapidly than any 
other body and will separate itself from the remainder of the swarm. Numer- 
ical models by Greenberg, et al. (1978) were the first to find this "runaway" 
path to planetary growth; subsequent calculations have confirmed these 
models (Wetherih & Stewart 1989, 1993; Ohtsuki & Ida 1990; Aarseth, 
Lin, & Palmer 1993). The basic mechanism behind runaway growth was 
elucidated by Wetherill & Stewart (1989) & (1993) Planetesimal velocities 
are determined by a balance among 1) stirring by gravitational scattering, 
2) stirring by inelastic collisions, 3) damping due to energy dissipation in 
inelastic collisions, 4) damping due to gas drag, and 5) energy transfer from 
large to small bodies via dynamical friction (Lissauer & Stewart 1993). The 
most critical factor is dynamical friction, which tends to lead to an equipar- 
tition of kinetic energy among the bodies (Binney & Tremaine 1987), with 
the largest bodies having the smallest random velocities and thus the largest 
cross-sections. Figure 6 illustrates the runaway process in which the largest 
body grows the fastest as long as the velocity dispersion remains low. 

Runaway growth can proceed until all the available material within the 
"accretion zone" (sometimes called the "feeding zone") of the planet is 
consumed (Wetherill 1980; Lissauer 1987, 1993). The accretion zone is an 
annulus over which the planet can exert its gravitational influence to per- 
turb nearby planetesimals into crossing orbits. The most common measure 
of the gravitational range of a planet of mass rrip a distance r away from a 
primary body of mass M^, is the Hill sphere (or Roche lobe radius) defined 
to be 

The origin of the cube root mass dependence comes from finding the region 
over which the gravity of the planet, Gm^/r"^, dominates the tidal force 
of the primary object, (GM*/r^)(rH/r). The radial size of the accretion 
zone has been estimated to be a numerical factor S ~ 4 times larger than 
the Hill sphere of the planet (Lissauer 1993). If the planet accretes all the 
available solid mass (with surface density Eg) in an accretion zone having 
width B th on either side, its final "isolation" mass will be 

mp = 27rr • 2SrH • Sg = A'KBr'^T.^ ( ^^ 
p MS y^j^^ 
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Figure 6. The evolution of the mass distribution of a planetesimal swarm between 0.99 
and 1.01 AU illustrating the rapid runaway growth of the largest bodies. The most mas- 
sive body becomes detached from the remainder of the swarm in less than 5 x 10'* years. 
See Wetherill & Stewart (1989). 






(27) 



24 



STEVEN P. RUDEN 



The radial spacing of planets that have each become isolated is ~ 2i?rH- 
For the MMSN, we find an isolation mass of 0.05M® at 1 AU and 1.4M© 
at Jupiter's distance. 
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Figure 7. Accretional evolution of a planetesimal swarm near 1 AU. The number of 
planetesimals are plotted per logarithmic interval in mass as a function of distance from 
the Sun. The evolution is clearly more rapid at 0.5 AU than at 1.5 AU, but the final 
mass of the protoplanets (~ 10^^ g) is not dependent on distance from the Sun. Run- 
away growth is evident as the largest bodies separate themselves from the rest of the 
planetesimal population. After 10® years there are more bodies with lower mass than in 
the terrestrial region of the Solar System. From Weidenschilling, et al. (1997). 



A recent detailed simulation of the accumulation of planetesimals in 
the terrestrial planet region has been performed by Weidenschilling, et 
al. (1997) that includes the processes mentioned above that control the 
velocity dispersion plus orbital perturbations caused by distant planetary 
bodies. The results of a simulation of the accumulation of 2 M® of ini- 
tially 4.8 X 10^® g (15 km) planetesimals spread between 0.5 and 1.5 AU 
is illustrated in Figure 7. These authors verified the critical role dynamical 
friction plays in producing rapid runaway growth. The first panel in Fig- 
ure 7 (t = 10^ years) clearly shows the evolution to larger masses occurs 



The Formation of Planets 25 

more rapidly at smaller radii where the Kepler frequency is larger. The 
other panels show that more massive bodies runaway and detach them- 
selves from the distribution of lower mass planetesimals. Runaway growth 
produces bodies with masses ~ 10^^ g in ;$10^ years, but the growth rate 
slows as the velocity dispersion of the population of lower mass bodies in- 
creases, and the runaway stalls. By 10^ years, about a dozen bodies with 
masses ~ 10^^ g have been produced, but they are lower in mass and more 
closely spaced than the terrestrial planets. The masses (and separations) 
found in the simulation are about a factor of three larger than the analytic 
isolation mass in equation (27). 

4.3. FINAL ACCUMULATION STAGES 

The outcome of runaway growth is self-limiting for two reasons. First, if 
the mass of the largest bodies becomes ~ 100 times the median mass of 
the continuum bodies, the velocity dispersion of the lower mass planetesi- 
mals is pumped up by gravitational scattering, and the runaway slows (Ida 
& Makino 1993; Weidenschilling, et al. 1997). Second, even if the velocity 
dispersion remains small, all the mass within the accretion zone will be con- 
sumed by the planet, and it will become dynamically isolated with a mass 
given by equation (27). The evolution following runaway is much longer 
term as gravitational encounters slowly perturb the bodies into crossing 
orbits and violent impacts occur. The models described in §4.2, which start 
from MMSN-type conditions, do not yield masses and spacings similar to 
those of the terrestrial planets, hence, longer term collisional evolution is 
necessary to produce a final model that is similar to the Solar System. Fig- 
ure 8 shows the results of six accumulation calculations by Wetherill (1988) 
that started with 500 bodies each of mass 2 x 10^^ g. The outcome of each 
simulation is a small number of planets after about 10® years. The differ- 
ences in the calculations illustrate the stochastic nature of the collisional 
accumulation. The long accumulation timescales do not present any serious 
difficulties for the terrestrial planets, which as noted in §1.1 could have com- 
pleted their evolution in a gas- free environment. However the giant planet 
cores had to form in less than the lifetime of the gas disk, which seems to 
require a rapid runaway growth process occurring in a disk with solid mass 
larger than the MMSN (Lissauer 1987, 1993; eq. [27]). 

4.4. FUTURE RESEARCH QUESTIONS 

— Can the Earth and the cores of giant planets be built directly from the 
runaway or is further long-term collisional growth needed? 

• How do you build the giant planet cores before the gas is dis- 
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persed? 

How does the final planetary system depend on the initial disk mass 
distribution? 

• Do higher mass disks produce more planets or fewer, but higher- 
mass, planets? 

Can rocky planets with masses comparable to Jupiter be built in pro- 
toplanetary disks? 



5. The Formation of Gas Giant Planets 

5.1. THE CORE-INSTABILITY SCENARIO 

The current model for the formation of gas giant planets is that they 
follow the path outlined above and begin as rocky cores. Mizuno (1980) 
showed that a gas envelope gravitationally bound to a rocky core can- 
not remain in hydrostatic equilibrium if the core mass exceeds a critical 
value ~ 10 — 15 M^, comparable to the core masses inferred for the giant 
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planets (Hubbard & Marley 1989; Chabrier, et al. 1992). Subsequent time- 
dependent calculations (Bodenheimer & Pollack 1986; Pollack, et al. 1996; 
see top panel in Fig. 9) have shown the nebular gas is rapidly accreted by 
the rocky core once its mass is above critical, leading to the conversion 
of a rocky planet into a gas giant planet. Before the critical core mass is 
reached, the radiated luminosity of the planet is supplied primarily by the 
gravitational energy of the accreted planetesimals. As the gravitationally 
bound gas envelope becomes more massive, planetesimal accretion is un- 
able to supply the luminosity needed to maintain the extended envelope and 
rapid (but not hydrodynamic) collapse of the envelope takes place, with the 
gravitational contraction of the envelope supplying the luminosity. As the 
envelope contracts, surrounding nebular gas flows in, and the gas accretion 
rate rises to become orders of magnitude larger than the solid accretion rate. 
The evolutionary timescale shortens; the gas envelope of Jupiter is accreted 
in ^10^ years (Pollack, et al. 1996). Detailed calculations show the value 
of the critical core mass {^ 10 — 30 M^) is only moderately dependent on 
orbital radius and nebular properties and depends most strongly on the 
planetesimal accretion rate (see Stevenson [1982] for a simple analytical 
calculation of the critical core mass) . 

In the core-instability model, giant planets are composite bodies with 
the rocky cores forming first and the gas envelopes accreting afterwards. 
If this model is correct, there is a logical simplicity to the formation of 
planetary systems - all planets are born by the accumulation of solids; 
giant planets are higher mass cores that happen to accrete large gas en- 
velopes. The core-instability model is not without problems, however. First, 
it is not obvious that ^lOM^ solid cores can form before the nebular gas 
is dispersed. Runaway growth to the critical mass requires solid densities 
about four times larger than in the MMSN (see eq. [27]). Forming the cores 
in ~ 10^ years requires a large net gravitational focusing factor (0^10^), 
which requires a very cold population of accreting planetesimals. Even if 
such cores can be formed at the orbit of Jupiter, it is not clear they can 
be formed so quickly at the orbits of Uranus or Neptune where the Kepler 
clock runs more slowly. Second, this model does not explain what halts 
the accretion of gas onto the planets. While it is possible the nebular gas 
was dispersed rapidly during the giant planet gas accretion phase, severe 
timing constraints (dispersal in ;$10^ years) make this unlikely. The most 
plausible explanation is that planetary gravitational tides become strong 
enough to "repel" nearby gas once the planet is sufficiently massive (Lin 
& Papaloizou 1993; see the next section). Finally, it is not clear from the 
internal structure of the giant planets in our Solar System that all four 
actually went critical. The large gas masses of Jupiter and Saturn seem to 
require their cores to have reached critical mass, but the smaller envelope 
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Figure 9. Core-Instability Model. Top panel from Bodenheimer & Pollack (1986) shows 
the core mass (solid) and envelope mass (dotted) as a function of time. Cases 6, 1 & 5, and 
7 have constant planetesimal accretion rates of 10~^, 10~^, and 10~^ M^/yr, respectively. 
Case 5 has grain opacity reduced by 50 to reflect the decreased dust-to-gas ratio after 
solids have sedimented from the gas. The envelope mass is very small until the critical 
core mass is reached (defined to occur when Mcore = Afonv). At that point rapid gas 
accretion ensues. Bottom panel from Pollack, et al. (1996) shows the core mass (solid), 
envelope mass (dot-dashed) and total mass (dotted) as a function of time. Rather than 
assuming constant planetesimal accretion rates (as in the top panel), this calculation uses 
more detailed and self-consistent planetesimal accretion rates. The initial core forms by 
runaway in less than l(f years after which the growth rate slows. The planet becomes 
critical after ~ 8 x 10 years, when rapid gas accretion begins. 



masses of Uranus and Neptune (Mcnv~0.5Mcorc) imply these planets did 
not. Shu, Johnstone, & Hollenbach (1993) have discussed how photoevapo- 
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ration of the disk by the ultraviolet flux from the central star preferentially 
removes disk gas beyond ^ 10 AU (the orbital distance of Saturn), which 
may explain the smaller gas mass reservoir available for the outer planets. 
Recent calculations of the core-instability model that adopt more realistic 
planetesimal accretion rates (Pollack, et al. 1996) indicate the planets ex- 
perience a long phase (-^ 7 x 10^ year) after core runaway in which the 
envelope mass is smaller than the core mass (see bottom panel in Fig. 9). If 
the nebula were removed during this long phase of "stalled" gas accretion, 
the resulting planets would be similar to the outer giants. 

5.2. TIDAL INTERACTION 

Gravitational interactions between orbiting mass points and disks has been 
considered by many authors (Goldreich & Tremaine 1979; Lin &: Papaloizou 
1979; Ward &: Hourigan 1989; Takeuchi, Miyama, & Lin 1996; see review 
by Lin & Papaloizou 1993). The essential physics can be elucidated by us- 
ing the impulse approximation (Lin & Papaloizou 1979). As gas flows by 
the planet interior (exterior) to its orbital radius, the planet's gravitational 
field causes a slight deflection in the gas streamlines causing the angular 
momentum of the gas to decrease (increase). The gravitational tidal in- 
teraction allows angular momentum to be transferred from the gas on the 
inside to the planet and from the planet to the gas on the outside. Be- 
cause orbital angular momentum increases outward in dynamically stable 
disks (the Rayleigh criterion; Drazin & Reid 1981), gas interior (exterior) 
to the planet that loses (gains) angular momentum in the tidal exchange 
must move farther inward (outward). In effect, the gravitational tides of 
the planet act to repel nearby gas. In a fluid dynamical context (Goldreich 
& Tremaine 1979), the tides of the planet excite spiral density waves that 
propagate away from the planet (in the absence of disk self-gravity). The 
waves are excited at Lindblad resonances, locations in the disk where the 
natural frequency of radial oscillation (the epicyclic frequency) is an integer 
multiple of the frequency at which the gas is forced by the planet.^ Spiral 
density waves carry angular momentum, which is deposited in the gas as 
the waves are damped by viscous stresses. Figure 10 illustrates the situation 
schematically. The total torque on the gas is found by summing the contri- 
butions from all the Lindblad resonances; it is negative (positive) for the gas 
interior (exterior) to the planet. The magnitude of the torque at either the 



^For Kepler disks, tlie Lindblad resonances are at radii satisfying ±m{Q{r) — Qp) = 
r2(r), where f2p is the orbital frequency of the planet and m is a positive integer. The 
positive (negative) sign is for inner (outer) resonances. 
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inner or outer resonances. 



Ttidal 



Ttidal, is approximately (Lin 

3 /^ \ 2 



f^K4 



rup 



Papaloizou 1993) 
(28) 



where ilp is the rotation frequency evaluated at the planetary radius rp, H 
is the vertical height of the disk and / ~ 0.2 is a numerical constant. This 
torque is deposited in the vicinity of the planet unless the viscosity is too 
small to damp the waves locally. The net tidal torque, which is the sum 
of the two nearly equal and opposite inner and and outer contributions, is 
difficult to calculate precisely but is roughly a factor H/r smaller than (28) 
(Ward 1986). 

GAP FORMATION IN PROTOPLANETARY DISKS 
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Figure 10. Formation of gaps in disks from Takeuchi, Miyama, & Lin (1996). Top 
Panel. Schematic illustration shows that gaps are formed from a balance between tidal 
torques that repel gas from the planet and viscous stresses that diffuse gas back into 
the gap. The tidal torque is actually deposited at the Lindblad resonances (indicated by 
arrows) . Bottom Panel. Schematic illustration of wave propagation in disks. Waves are 
excited at the inner and outer Lindblad resonances and propagate away from the planet. 
As they damp, the angular momentum they carry is deposited into the gas as a tidal 
torque. A net outward flux of angular momentum results. 

The action of the planetary tidal torques, which tend to push gas away 
from the planet, is opposed by the viscous stress of the gas, which tends 
to diffuse gas back toward the planet. The viscous torque on the disk gas, 
Tvis, is (Lynden-Bell & Pringle 1974) 



37rSz^ Q.r 
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-)■ 

r I 



(29) 
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using equations (3) and (7) to get the rightmost form. The criterion for the 
planet to open an annular gap at its orbital radius is that the gravitational 
tidal torque exceed the viscous torque, Ttidai~rvis. Using equations (28) 
and (29) the minimum planet mass, rrip, needed to open a gap is (Takeuchi, 
Miyama, & Lin 1996) 




(30) 

For typical solar nebular parameters this gives a mass of ~ 75M^, slightly 
less than the mass of Saturn. It thus seems reasonable that Jupiter and Sat- 
urn (but not Uranus or Neptune unless the disk was extremely cold) were 
able to form gaps in the surface density distribution, reducing the local 
gas density to small values and possibly terminating the rapid gas accre- 
tion phase. I note however that some recent hydrodynamical simulations of 
the gap formation process though have questioned whether gap formation 
completely halts gas accretion onto the planet. The simulations indicate 
that even if a low density gap is formed, the planet remains connected to 
the gap walls by arc-like streams through which gas can flow (Artymowicz 
& Lubow 1996, personal communication 1998). Clearly far more work is 
needed to determine whether the formation of a gap completely halts gas 
accretion, and if not, what determines the final mass of the giant planets. 

5.3. TIDAL MIGRATION 

An essential consequence of the tidal exchange of angular momentum be- 
tween the planet and disk is that a net torque can act on the planet, causing 
its orbital radius to change with time. This process is known as orbital or 
tidal migration (Lin & Papaloizou 1986; Ward 1997). The angular momen- 
tum of the planet (on a circular orbit) is mp^pVp = m,p(GM*rp)^/^ and 
writing the net tidal torque exerted by the planet on the gas disk as T^et , 
the orbital evolution of the planet is governed by 

d_ 
di 
or equivalently 



(^mpQprlj = -Tnet, (31) 

Tnet ■ (32) 



drp _ 2 

dt TTipilprp 



The orbital migration timescale is 
t 



t-rnig 




— , (no gap) (33) 
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where I have used equation (28) reduced by H/r to estimate the net torque 
(Ward 1986). This equation is appropriate for lower mass planets that do 
not tidally truncate the disk and form gaps. For planets massive enough to 
form gaps, the surface density in their vicinity is reduced dramatically, and 
equation (33) does not apply. Numerical calculations of the combined gap 
formation and tidal migration process (Lin & Papaloizou 1993) show that 
the inner and outer torques nearly cancel and the evolution time becomes 
the local viscous diffusion time, equation (8), applied at the orbit of the 
planet. Figure 11 illustrates a self-consistent calculation of gap formation 
and tidal migration in a protoplanetary disk. 




Figure 11, Tidal migration of a protoplanet in an evolving nebula from Lin & Pa- 
paloizou (1993). Surface density, S, is plotted versus distance from star, r, as a function 
of time, t. All variables are scaled and dimensionless. Notice the surface density de- 
creases with time (due to accretion) , and the outer disk radius increases. The position of 
the planet is marked by a dot. The mass of the planet is large enough that tides over- 
come viscous stresses, and a gap forms near the planet (the region where S — > 0). In 
this example, the planet is near the outer disk radius and migrates outward as it absorbs 
angular momentum from the inner disk. Eventually the inner disk will be accreted by the 
star, and the planet will migrate inward. The migration timescale is the viscous timescale 
of the disk. 



The tidal migration timescale is surprisingly short. Using values ap- 
propriate to the MMSN, the timescale (33), appropriate for planets not 
forming gaps, is tmig ~ {Mj, / m-p){r / Pi\]) years, which is about 3 x 10^ years 
for the Earth. For planets massive enough to form gaps, the ratio of the 
viscous time at their orbital distance to the overall viscous time of the disk 
(having outer radius Rd) is ~ {r^/ RdY'"^ , from equation (8) assuming H/r 
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is roughly constant. Because Jupiter and Saturn were well within the outer 
radius of the primordial solar nebula (rp <C Rb), their migration times were 
much shorter than the ~ 10^ year lifetime of the disk as a whole. Figure 12 
shows a calculation of the orbital migration of Jupiter and Saturn in which 
a migration time of ~ 10^ years is found. 

5.4. THE SURVIVAL OF GIANT PLANETS 

Tidal migration offers a neat way to explain the recently discovered Jupiter- 
mass planets having orbital distances of ^0.1 AU from their stars (Marcy & 
Butler 1998; also see the review by Marcy in this volume) if a mechanism 
can be found to halt the migration before the planets are swallowed by 
the star. The key is to find an outward torque that can counterbalance the 
inward tidal torque. Lin, Bodenheimer & Richardson (1996) have suggested 
that as the planet nears the central star, it will raise tidal bulges in the 
star that will transfer angular momentum from the spin of the star to 
the planetary orbit. This spin-orbit coupling can be effective at keeping 
the planet beyond the corotation radius (the distance in the disk where 
the Kepler orbital frequency equals the stellar rotation frequency). These 
authors also suggest that a star with a sufficiently strong magnetic field will 
have a magnetosphere that truncates the disk at a distance of several stellar 
radii. Once the planet migrates inside the inner disk edge, the inward tidal 
torques vanish because the surface density is negligible. Either mechanism 
can produce planets in orbits near the star. 

Trilling, et al. (1998) have considered a mechanism involving mass loss 
from the planet. The size of the planetary Roche lobe (c./. eq. [26]) is pro- 
portional to the distance from the star. As a giant planet migrates inward, 
it can fill its Roche lobe, leading to a transfer of mass from the planet to 
the star. Conservation of angular momentum of the system requires the 
planet to move outward. Stable mass transfer is achieved when the planet 
orbits at a distance where its radius equals its Roche radius. Planets in 
small orbits near stars can be produced if the disk is dissipated before the 
planet loses all its mass. 

Our own Solar System with planets distributed from 0.4 - 30 AU does 
not shown signs of appreciable orbital migration, however. The question in 
our Solar System is to understand how the large orbital migrations appar- 
ently observed in other planetary systems were avoided. Because the mi- 
gration timescale becomes the viscous timescale for planets massive enough 
to form gaps, a giant planet can avoid appreciable migration if it is born in 
a low viscosity nebula having evolutionary times comparable to or longer 
than the disk dispersal time. Before the planet has time to migrate far, the 
nebular gas is dispersed, and the migration is halted. Denoting the disk 
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Figure 12. Tidal migration in the solar nebula models of Ruden & Lin (1986). Surface 
density in g/cm^ is plotted versus distance from the Sun in AU. Jupiter is assumed to form 
instantaneously with its present mass and orbital radius at a time in the nebular evolution 
when it is just massive enough to form a gap (eq. [30]). Saturn is assumed to form coevally. 
Top panel. The solid curve is 10'^ years after Jupiter formation showing the gap (the 
region where E drops to zero) is formed rapidly. The dashed curve is 4 x 10* years later, 
showing Jupiter has migrated inward to 1 AU. Bottom panel. The solid curve is also 
10"^ years after Saturn formation. Saturn does not form a gap, but the surface density is 
reduced there by 30%. The dashed curve is 1.5 x 10^ years later. At this later time, the 
surface density in the gap is not zero but has been reduced by a factor of 5. The average 
inward migration speeds are about 40 cm/s for either planet. 



dispersal time by tdispi we can estimate from equation (8) that this requires 
the viscous alpha parameter to be ^(r/if)^/(r2tdisp)- Using r /H ^ 25 and 
idisp ~ 10 years, I find q;;S10~^ is sufficiently small that migration of 



Jupiter will not be severe. It is not unreasonable that giant planets form in 
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a low viscosity environment. The lower panel in Figure 9 shows a model for 
the formation of Jupiter in which rapid gas accretion does not begin until 
t ~ 8 X 10^ years. In this model, the newly formed Jupiter would find itself 
in a nebular environment that had undergone significant evolution. The 
disk mass at that epoch could be quite small, with the corresponding low 
surface densities leading to small tidal torques and low migration rates. In 
fact, given that the MMSN contains only 10 - 20 Jupiter masses, formation 
of the giant planets at a late epoch when the gas mass may be appreciably 
less than in the MMSN may require the giant planets to accrete a large 
fraction of the remaining disk mass. 

5.5. FUTURE RESEARCH QUESTIONS 

— Is the core-instability model the only way to make giant planets? 

• Did the dispersal of the disk prevent Uranus and Neptune from 
going critical? 

— What determines the final masses and orbital distances of the giant 
planets? 

• Does gap formation halt gas accretion onto giant planets? 

• What prevented significant orbital migration from occurring in 
our Solar System? 

— How does the planetary mass spectrum merge into the brown dwarf 
mass spectrum? 

• Do brown dwarfs form in a fundamentally different manner than 
planets, e.g., from collapse/fragmentation of the disk rather than 
from accumulation of dust? 
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